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Abstract The luminosity for a WASA-at-COSY experiment involving the pd reaction at 2.14 GeV proton- 
beam energy is determined by the forward pd elastic scattering, which yields an average beam-on-target value 
of [5.2±0.3(stat)±0.3(syst)] x 10 30 s _1 cm -2 . In addition, the forward pd elastic-scattering angular distribution 
is obtained with four-momentum transfer squared —t between 0.16 (GeV/c) 2 and 0.78 (GeV/c) 2 at this beam 
energy, which is compared with other experimental data and the pd double scattering model. 
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1 Introduction 

The WAS A detector facility 0, a nearly 4n mul- 
tidetector system equipped with an internal H 2 /D 2 
frozen-pellet target [2j , is now operated at the COoler 
SYnchrotron COSY- Jiilich [3| which delivers proton 
and deuteron beams with momenta up to 3.7 GeV/c. 
It serves for the investigation of hadronic processes 
and systems, such as light meson production, exotic 
hadronic states, symmetries and their breaking, fol- 
lowing pp, pd and dd reactions Q. WASA is well 
suited to detect both charged and neutral decays, 
e.g. of the light scalar mesons a //o(980) produced 
in such reactions. A first test measurement has been 
performed in order to obtain the yet unknown light- 
scalar production cross sections in pd — > 3 AX pro- 
cesses at a beam energy of T p = 2.14 GeV from the 
measurement of its strong decays. We describe here 
the luminosity determination for this experiment. 

The luminosity is defined as the number of beam 
particles passing through the target per unit time 
multiplied by the number of atoms in the target per 



unit area. In the case of a stored beam and an internal 
pellet target, the luminosity can be inferred from the 
number of beam particles, the beam-revolution fre- 
quency, the number of atoms in one pellet, and the 
beam-pellet overlap factor. The parameters of the 
COSY beam and the pellet target for this particular 
experiment are listed in Table [TJ 

The proton beam in the 184 m long storage ring, 
circulating with a revolution frequency / related to 
its energy T p , has an average beam intensity N c and 
a diameter 2R C at the interaction region. The num- 
ber of beam particles per unit time and area, called 
current density, is obtained from the formula: 

j = (N c f)/(nR%) . (1) 

On the other hand the number of atoms in one D 2 - 
pellet can be calculated from its radius R t and density 
p with: 



iV t (D) 



(2) 



where n = 2 atoms/molecule for Deuterium, M ss 4 
g/mol is its molar mass, and N A — 6.022 x 10 23 
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molecules/mol is the Avogadro constant. Inserting 
the numbers from Table [TJ the average luminosity for 
this experiment is estimated as I ~ 5 x 10 30 s _1 cm -2 
which should be regarded as an order-of-magnitude 
estimate, since the beam diameter and the pellet rate 
vary as a function of time. 

Table 1. Beam and target parameters. 



a much smaller width of a = 2.0° is obtained. 
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Parameters 


Values 


Units/remarks 


Beam intensity (Nq) 


1.7 X 10 9 


particles 


Revolution frequency (/) 


1.5564 MHz 


2.14 GeV 


Beam size (2Rc) 


6 mm 


3-6 mm 


Current density (J) 


9.3 X 10 15 


particles / (s • cm 2 


Pellet diameter (2R t ) 


35 fim 


25-35 /im 


D2 density (p) 


0.162 g/cm 3 


triple point 


Atoms/pellet (JVt(D)) 


1.10 x 10 15 


D 2 


Pellet frequency 


6 kHz 


5-12 kHz 


Pellet velocity 


70 m/s 


60-80 m/s 


Pellets in beam 


0.5 


overlap factor 



2 Identification of pd elastic scattering 

Besides the estimate outlined above, it is useful 
for luminosity determination to record a monitoring 
reaction during the experiment, such as elastic scat- 
tering for which the corresponding cross sections are 
well known in dependence of both the beam energy 
and the scattering angle. In the data acquisition 
(DAQ) the trigger settings for two charged tracks 
within WASA were applied to select elastic-scattering 
events. Both pd elastic scattering and pp quasi-elastic 
scattering, with the neutron being a spectator parti- 
cle, could be recorded by the same trigger condition. 

For the scattering processes the data analysis 
starts with the selection of one charged track in the 
forward detector (FD, covering emission angles of 
9 = 3—18°) and the second one in the central detector 
(CD, 9 = 20- 170°). After a FD-CD time correla- 
tion cut, the coplanarity of the two charged tracks is 
exploited. The corresponding azimuthal-angle corre- 
lation is plotted in Fig. Q] For the azimuthal-angle 
correlated events, the polar-angle correlation of the 
charged FD and CD tracks is also depicted in Fig. [T] 
It is seen that most events are located around the pp 
elastic line (determined from a MC simulation) and 
thus stem from pp quasi-elastic scattering. The MC 
simulations are carried out with the Pluto event gen- 
erator [6| and the WASA detector simulation pack- 
age based on GEANT3. The peak around 180° of the 
azimuthal-angle correlation has a width of a = 6.9° 
in data, which is reproduced by the simulation of pp 
quasi-elastic scattering, while for pp elastic scattering 
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Fig. 1. Angular correlation of two charged 
tracks for the azimuthal (upper panel) and po- 
lar (lower) angles. The dashed line presents 
the MC simulation of pp quasi-elastic scatter- 
ing from pd reaction. Line A depicts the kine- 
matics of pp elastic scattering; Lines B and 
C are for pd elastic scattering with the proton 
emitted in forward and backward direction, re- 
spectively. 

The pd elastic-scattering events, with the proton 
detected in the FD and the deuteron in the CD, can 
be identified on top of a huge background from the 
above mentioned pp quasi-elastic scattering. For the 
charged track in the CD, the AE-vs. -momentum dis- 
tribution is presented in Fig. [2] together with the ex- 
pected proton and deuteron lines obtained from a MC 
simulation. A weak cut is applied, which then allows 
one to extract the pd elastic-scattering events from 
the strongly reduced pp quasi-elastic background. For 
the deuteron candidates the agreement between the 
polar-angle value measured directly in the CD and 
the one calculated from the other charged track in the 
FD is checked on an event- by- event basis. Figure [2] 
shows that a clear peak of forward pd elastic scatter- 
ing appears around zero angular difference. We did 
not further analyze the backward pd elastic scatter- 
ing with the deuteron in the FD and the proton in 
the CD since there are very few events around Line 
C shown in Fig. [1] 
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Fig. 2. Upper panel: Energy deposited in the 
plastic scintillator vs. momentum recon- 
structed in the solenoid field of WASA. Lower 
panel: Difference of the measured polar an- 
gle of the charged track in the CD and the 
one calculated from the measured polar angle 
in the FD assuming forward pd elastic kine- 
matics. The background under the pd elastic- 
scattering peak has been fitted by a fourth or- 
der polynomial and is shown by the green line 
(color online). The remaining signal is also 
indicated in the figure by a Gaussian distribu- 
tion. 



3 Normalization of pd angular distri- 
bution 



0.78 (GeV/c) 2 . The number of pd elastic-scattering 
events is extracted from the number of total events in 
the gaussian ±3cr region with the relative area ratio 
of the gaussian peak to the polynomial background, 
shown in Fig. [2 and the statistical error is taken from 
the number of total events and the fitting error of 
the polynomial background. The geometric accep- 
tance, reconstruction and cut efficiencies have been 
determined as a whole from the MC simulation for 
each polar-angle interval as listed in Table [2] In ad- 
dition, the trigger pre-scaling factor is 2 000 for the 
first 35 runs and then 4 000 for the later 123 runs. 
The DAQ lifetime correction is estimated with a rel- 
ative error 12.4%. Both, the statistical error and the 
DAQ lifetime-correction error, are listed in Table [2j 



Table 2. Forward pd elastic-scattering angular 
distribution at 2.14 GeV from our data. The 
last column (dN/dfi)i a b includes the statisti- 
cal error of the pd elastic-scattering events and 
the error of the DAQ lifetime correction. 



proton, lab 


-t 


efficiency 


(dN/dQ) lilh 


C) 


(GeV/c) 2 


(96) 


(xl0 6 /sr) 


8.5 


0.186 


10.4 


5640±190±699 


9.5 


0.231 


42.0 


884±31 ±110 


10.5 


0.281 


49.3 


532±24 ±66 


11.5 


0.336 


51.1 


314±16 ±39 


12.5 


0.394 


51.2 


238±12 ±30 


13.5 


0.457 


51.6 


175±9 ±22 


14.5 


0.523 


52.2 


210±11 ±26 


15.5 


0.594 


52.4 


139±8 ±17 


16.5 


0.668 


50.2 


138±8 ±17 


17.5 


0.745 


38.2 


118±10 ±15 



The normalization of pd elastic-scattering angular- 
distribution (dN/dft) lah to the differential cross sec- 
tions (da/dfl)^ is described by the formula: 



(3) 



where L int is the integral luminosity for the corre- 
sponding beam time. 

The forward pd elastic-scattering angular distri- 
bution from our data with four-momentum trans- 
fer squared —t, given in Table [2j is obtained in a 
polar-angle interval of 1° in the FD with the pd 
elastic-scattering peak observed starting from 8° to 
18° in the lab frame. The —t is the Lorentz-invariant 
momentum transfer ranging from 0.16 (GeV/c) 2 to 



The forward pd elastic-scattering differential cross 
sections from 2.0 GeV data ll| with four-momentum 
transfer squared —t between 0.35 (GeV/c) 2 and 
1.0 (GeV/c) 2 are fitted with the exponential func- 
tion a-e 6 ' - **, shown in Fig. El which has a slope value 
&i = -1.87±0.09, while a slope b 2 = -1.84±0.50 is ob- 
tained from fitting our data independently. When the 
slope value is fixed as b = —1.87 for both cases, the fit- 
ting procedure gives the coefficient = 710±11 fib/ sr 
for 2.0 GeV data and a 2 = 465±26 /ib/sr-pb -1 for our 
data, respectively. Then, a scaling factor is applied 
to normalize our data to the same height of 2.0 GeV 
data, which leads to the integral luminosity. 
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Fig. 3. Normalization of the pd elastic- 
scattering angular distribution from our data 
to the differential cross sections from 2.0 GeV 
data with —t between 0.35 (GeV/c) 2 and 



1.0 (GeV/c) 2 . 
nential. 



The fit function is an expo- 



The pd elastic-scattering differential cross sec- 
tions for protons in the forward direction with four- 
momentum transfer —t less than 1.0 (GeV/c) 2 are 
plotted in Fig. U for proton-beam energies ranging 
from 0.425 GeV to 11.90 GeV BO- 
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• 425 MeV [7] 

• 582 MeV [8] 
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Fig. 4. Compilation of forward pd elastic- 
scattering differential cross sections vs. four- 
momentum transfer —t including our data af- 
ter the normalization. 

The forward pd elastic-scattering angular distri- 
bution can be explained by nucleon-nucleon single 
and double scattering and their interference for the 
deuteron as a double scatterer [3) El]- A further 
theoretical investigation, which takes the d-state ad- 
mixture fully into account fl3| , gives a better descrip- 
tion of the shoulder near —t ~ 0.44 (GeV/c) 2 . Our 
data points are in good agreement with the model 
calculations around this sensitive region (see Fig. [5]) 
which provides further confidence in the normaliza- 
tion procedure. 
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Fig. 5. Comparison of the forward pd elastic- 
scattering angular distribution available from 
our data with the model calculation [l^] after 
the normalization. 



4 Result and error estimate 



The average luminosity in the units of (s~ 1 cm~ 2 ) 
can be obtained from the integral luminosity divided 
by the effective beam time, which discards the rest 
time between the cycles. There are 274 runs in 
this experiment with a data-taking time of about 
80 hours, of which the first 158 runs have the valid 
trigger condition for the forward pd elastic scattering 
with the effective beam time 1.19 x 10 s seconds. The 
average luminosity is obtained from this part of the 
beam time. 

In the normalization of the angular distribution, 
the beam energy difference between our data and 
2.0 GeV data is not considered, which could intro- 
duce a small correction. The pd elastic-scattering 
differential cross sections in the lab frame with four- 
momentum transfer squared —t at 0.44 (GeV/c) 2 are 
plotted in Fig. [6] for proton-beam energies ranging 



im encrgic 

from 0.425 GeV to 11.90 GeV @, H Q [l2|. Af- 
ter comparing the agreement between data and three 
kinds of fit functions, the quadratic one is used to esti- 
mate the correction factor with a value c = 1.06±0.07 
for the differential cross section when the beam en- 
ergy increasing from 2.0 GeV to 2.14 GeV. 
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Fig. 6. Fit of pd differential cross sections at 
the same —t value vs. proton-beam energy. 
Three kinds of fit functions are compared. 

Then, the integral luminosity is deduced as 

L int (158 runs) = — ^~ 
c-ai 

= 0.62±0.03(stat)±0.04(syst) pb" 1 , (4) 

with the statistical error from a 2 and the systematic 
error from c and a 1 . This yields the average luminos- 
ity during beam-on-target times as 

L = [5.2±0.3(stat)±0.3(syst)] xlO^'s^cnT 2 . (5) 

Finally, it is possible to estimate the total inte- 
gral luminosity with the help of another monitoring 
trigger defined as two charged tracks in the central 



detector, which was available during the whole beam 
time. These two monitoring triggers have almost the 
same rates after pre-scaling, and their ratio is more or 
less stable within a relative error of 6.5%. Including 
the above relative error 12.4% from the DAQ lifetime 
correction, the ratio of the total integral luminosity to 
the partial integral luminosity of Eq. [4] is 1.54 ±0.22, 
which yields 

L int (274 runs) = 
0.95±0.05(stat)±0.15(syst) pb _1 . (6) 

5 Summary 

The aim of this experiment is to measure the light 
scalar meson a // (980) production in pd^ 3 AX re- 
actions, while the pd elastic scattering is measured 
in parallel as a reference reaction to determine the 
luminosity. As the effective beam time is the same 
for both aim and reference reactions, the value of the 
integral luminosity in Eq. [5] is used for the evalua- 
tion of the a // (980) production cross sections in pd 
reactions. 

We gratefully thank the COSY operators, the tech- 
nical and administrative staff at the Forschungszen- 
trum Jiilich as well as all the members of the WASA- 
at-COSY collaboration for their support. 
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